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We show that the entanglement entropy and alpha entropies corresponding to spatial polygonal 
sets in (2 + 1) dimensions contain a term which scales logarithmically with the cutoflt. Its coefficient 
is a universal quantity consisting in a sum of contributions from the individual vertices. For a free 
scalar field this contribution is given by the trace anomaly in a three dimensional space with conical 
^ ' singularities located on the boundary of a plane angular sector. We find its analytic expression as a 

O , function of the angle. This is given in terms of the solution of a set of non linear ordinary differential 

' equations. For general free fields, we also find the small-angle limit of the logarithmic coefficient, which 

is related to the two dimensional entropic c-functions. The calculation involves a reduction to a two 
^is^ . dimensional problem, and as a byproduct, we obtain the trace of the Green function for a massive 

scalar field in a sphere where boundary conditions are specified on a segment of a great circle. This 
also gives the exact expression for the entropies for a scalar field in a two dimensional de Sitter space. 

> 
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1 Introduction 

Suppose we are interested in the physics of events locahzed in a region V of the space. The Hilbert space 
of states can be decomposed accordingly as a tensor product H = Hv ^ TC-v of the spaces of the states 
rH I localized in V and in the complementary region —V. Take now the vacuum l^*) as a global state of the 

system, with density matrix po = j^*) (^fj. The state pv relevant to the algebra of operators acting on 
Hv follows from the partial trace of po over the complementary Hilbert space H-v- This gives the local 
' reduced density matrix 

^ : Pv = trw_. If) {^\ . (1) 

The global state po is generally entangled in the bipartite system Hv ^ 'H-v and in consequence this 
matrix is mixed. The corresponding entropy 

S[V)^~iv{pv\ogpv) (2) 

is usually called entanglement or geometric entropy. 

The entanglement entropy is one of the most prominent candidates to explain the intriguing entropy 
of the black holes [1]. However, in this proposal the role of quantum gravity is fundamental to produce 
a finite entropy, and the whole subject is still controversial. On the other hand, this and other measures 
of entanglement have also been extensively studied in condensed matter and low dimensional systems, 
partially motivated by advances in quantum information theory and the density matrix renormalization 
group method. As a result it was uncovered that a variety of phenomena such as quantum phase transitions 
have an interesting correlate in the entanglement properties of fundamental states 2[. From the point of 
view of quantum field theory (QFT) the function S{V) can be considered as a non local variable with 
interesting non perturbative properties [3J, which can be defined for any theory disregarding the field 
content. In this context, we can mention among the applications the description of topological order [3], 
and the renormalization group irreversibility in two dimensions [51 [B] . 



*e-mail: casini@cab.cnea. gov. ar 
t e-mail: hucrta@cabtep2.cnea. gov. i 



1 



In the continuum limit described by a quantum field theory the entanglement entropy is divergent 
due to the presence of an unbounded number of local degrees of freedom. The divergent terms must be 
proportional to quantities which are local and extensive on the boundary of V. This can be seen as a 
consequence of the local nature of the ultraviolet divergences, and that the boundary is shared between V 
and —V which have the same entropy (for any pure global state) 

S{V)^Si~V). (3) 

On a technical level this characteristic of the divergent terms is due the fact that the entanglement entropy 
is the variation of the euclidean free energy with respect to conical singularities located on the boundary 
of y [7|. Thus, in d spatial dimensions, we expect to have an expansion of the form 

S{V) - gd-i[dV] e-^''-^^ + gd-2[dV] e-^'^-^) ^ + g^[dV] log(eA) + So{V) , (4) 

where Sq{V) is a finite part, e is a short distance cutoff, and the gi are local and extensive functions on 
the boundary dV, which are homogeneous of degree i (in a different context a similar expansion was used 
in [8]). The leading divergent term coefficient gd-i[dV] is proportional to the d—1 power of the size of V. 
This was noted since the earliest papers on the subject [H [9] and is usually referred to as the area law for 
the entanglement entropy. However, strictly speaking, gd-i depend on the regularization procedure and 
it is not proportional to the area if this later is not rotational invariant. For example, if we draw a square 
of side L on a two dimensional lattice, we would have gi{dV) ~ Ci-^, where e is the lattice spacing, but 
where the dimensionless constant Ci strongly depends on the relative angle between the square sides and 
the lattice symmetry axes. Moreover, all the terms gi with i > 1 are not physical within QFT since they 
are not related to continuum quantities. 

On the contrary, the dimensionless coefficient go{V) of the logarithmic term is expected to be universal. 
Logarithmic divergent terms in the entropy have been previously found in four dimensional black hole 
space-times [lOj . They are present generically in even dimensions for sets with smooth curved boundaries. 
This follows from the heat kernel expansion for conical manifolds with smooth singularity surface [TT] . 

In this work we show that there is also a logarithmic term in three dimensions for sets V with non- 
smooth boundary. In particular, we consider the case of spatial polygonal sets (in a different scenario, 
a logarithmic contribution to the entanglement entropy was also reported in Since go{V) is di- 

mensionless, extensive and local on the boundary we conclude that for V a polygon it must be of the 
form 

5o(V^) = E^(^')' (5) 

where the sum is over all vertices Vi and Xi is the vertex angle. On general grounds one also expects 
point-like vertex induced logarithmic terms in any dimensions. 

All these considerations apply to the alpha-entropies (or Renyi entropies) as well 

So.{V) - -^log(trp^), S{V) = lim Sa,{V) . (6) 

L — a a^l 

These measures of information behave very much like the entanglement entropy, as discussed in several 
papers [T31 [T3] . For integer a = n they are also more suitable for explicit calculation, since the traces 
trpy involved can be represented by a functional integral on an n-sheeted space with conical singularities 
located at the boundary of [J. The entanglement entropy follows from Sn by analytical continuation. 
Thus, in 2+1 dimensional theories the alpha-entropies contain a logarithmic term 

'^"liog = IIs"(2^01og(eA), (7) 

analogous to the one in (|3|) for the entropy. We also have s{x) — lima^i Sa{x). In eq. ([7]) A is a parameter 
with the dimensions of an energy, depending on V and on the particular theory. For a massless field it is 
the inverse of any typical dimension R ol V and when the mass dominates, MR 3> 1, it can be taken as 
A = M. We find the analytic expression of Sn{x) for a free scalar field, and check our results by numerical 
simulations on a two dimensional lattice. 

The outline of the paper is the following. In Section 2, we obtain the analytic expression for trpy for 
the local density matrix associated to a plane angular sector for a massive scalar field in 2 -|- 1 dimensions. 
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The calculation is reduced to the one of the trace of the Green function on a two dimensional sphere, where 
boundary conditions are specified on a segment of a great circle. We use the method introduced in [M] to 
obtain this trace exactly. Finally, we identify the logarithmic coefficients s„(x) of ([7]). We also compare 
them with the results given by simulations in a square lattice finding a perfect accord. In Section 3, we 
find the small angle limit of s{x) for free fields and relate it to the two dimensional entropic c-functions. 
Finally, in Section 4 we present our conclusions. We have also included in the Appendix A a detailed 
derivation of the formulas concerning the Green function in the sphere with a cut presented in Section 
2.1, and in Appendix B, the method we have used to compute numerically the geometric entropy in a two 
dimensional lattice. 



2 Vertex induced logarithmic terms in the entropy 

To compute the coefficients s„ (x) of the logarithmic term in ([7|) for a scalar field we choose for convenience 
to study the entropy associated to a set V given by a plane angular sector of angle x. 

The traces tipy involved in ([5]), with a = n G Z, can be represented by a functional integral on an 
n-sheeted three-dimensional Euclidean space with conical singularities located at the boundary of the set 
V [71 [13]. To be explicit, calling u and v to the upper and lower faces of the plane angular sector, the 
replicated space is obtained considering n copies of the three-space cut along the angular sector, and sewing 
together the upper side of the cut u*"' with the lower one v'''^^ , for the different copies k = l,...,n, and 
where the copy n + 1 coincides with the first one. The trace of p" is then given by the functional integral 
Z[n] for the field in this manifold, 

trp" = (8) 

Then, following [131 114| we consider a free massive complex scalar field on this manifold, and map 
the problem to an equivalent one in which we deal with n decoupled and multivalued free complex scalar 
fields. First we arrange the values of the field in the different copies in a single vector field living in a three 
dimensional space, 

/ 01 (^) 



$ = 



(9) 



where ^i{x) is the field on the I copy. Note that in this way the space is simply connected but the 
singularities at the boundaries of V are still there since the vector $ is not singled valued. In fact, crossing 
the plane angular sector from above (side u) or from below (side w), the field gets multiplied by a matrix 
T or T^^ respectively. Here 

/ 1 \ 

1 
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1 

/ 



(10) 



which has eigenvalue|l| e*"^'^, with fc = , (n — 1). Then, changing basis by a unitary transformation 
in the replica space, we can diagonalize T, and the problem is reduced to n decoupled fields 0fc living on a 
single three dimensional space. These fields are multivalued and defined on the euclidean three dimensional 
space with boundary conditions imposed on the two dimensional set V given by 



~^l{T) = e'^~^l(T) A: = 0, ...71-1, fey. 



(11) 



Here 0^ and i/)^, are the limits of the field as the variable approaches V from each of its two opposite sides 
in three dimensions. In this formulation we have 



log(trp^) = ^ogZ[kl7 



k=0 



(12) 



^ There is an important difference with the free fermion case studied in 1131 , where the eigenvalues of T are e* n ^'^ , 
k = — (n — l)/2 , (n — l)/2. This is because in the fermionic case it is T" = —1. 



with 



3 



where Z[a] is the partition function corresponding to a field which acquires a phase e*^'^" when the variable 
crosses V. We note that the coefficient of the logarithmic divergent term in log Z[k/n] corresponds to the 
integrated trace anomaly J dr^Tj^{r), or equivalently the 03 — 0^/2/ {'^t^)'^^'^ coefficient of the heat kernel 
expansion [17j . for a massless scalar in three dimensions with the boundary conditions (jlip . 

Our calculation of Z[a\ is based on the relation between the free energy and the Green function for a 
free massive scalar 

aM^logZ[a] = - j dr^Gaif, ^ ■ (13) 

Here Ga{r 1,7^2) is the Green function for a complex scalar of mass M in three dimensions subject to the 
boundary conditions (|lip . To be explicit, we have 

+Af2)Ga(fi,f2) - d{n~r2), (14) 
lim Ga(ri +er),r2) = e'^^'^ lim Ga(fi - £??, fa) , n e V , (15) 

where 77 is orthogonal to the plane of V. 

The Laplacian and the boundary conditions allow the separation of angular and radial equations in polar 
coordinates. Using standard methods we arrive at the expression (see for example the similar calculation 

mm) 

Ga{ri,r2)=Y. J dX ^^^^^^ M0i,^iWA02,V2) '^'^ , (16) 

where J is the standard Bessel function. Here the sum is over the normalized eigenvectors V'fC^i f) of the 
angular equation 

l^ni'u ^ + , (17) 

where Aq is the Laplacian on the sphere with domain given by the functions satisfiyng the boundary 
conditions inherited from (jlSp . Specifically we can choose the cut on the equatorial plane 

lim Vi^(7r/2 + e,ip)^ e'^^ra j-^^ ^..(71/2 - e,ip) , Lp £ [~x/2, x/2\ . (18) 

e— >0+ e— >0+ 

The functions are the spherical harmonics if a = 0, otherwise they are known to satisfy Lame differential 
equations [16j . Also, the values of u are not integer if a 7^ 0. The precise expressions for and v will not 
be relevant in what follows. 

Taking the trace / dr^Ga{r,r) in eq. (fTB)) gives 



logZ[a] ^ --^ ^(^ + 1/2) = -^tJ-An + ^ . (19) 



2M2 ' ' V 4 

Though this expression is divergent, the piece we are interested in, which is the one dependent on the 
angle x, is finite. 

To proceed, we find convenient to express the trace of the square root of the operator in in terms 
of the corresponding Green function. This is done by using the expression for the powers of an elliptic 
operator O in terms of the resolvent given in [18] , 

O'' = — I X<{0 -Xy^dX, (20) 

where F is a curve depending on the particular operator O. In the present case F begins at infinity, pass 
along the negative real axes on the upper complex plane, encircles the origin and goes back to infinity on 
the lower half of the complex plane along the negative real axes. This gives for the trace 
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2.1 Green function on a sphere with a cut 

The problem is then reduced to the calculation of the trace of the two dimensional Green function on a 
sphere with a cut of angle x, where the boundary conditions (HH]) are imposed. In a previous paper [13] we 
have solved the analogous problem on the plane with boundary conditions imposed on an interval, using a 
generalization of a method originally developed in ^19j . It essentially consists in exploiting the symmetries 
of the Helmholtz equation even in the presence of symmetry breaking boundary conditions by analyzing 
the behavior of the Green function at the singular points. 

Following the recipe of jT4] step by step (with the additional algebraic complications corresponding 
to the spherical case) we find the analytic expression for the trace of the Green function as a solution of 
a system of ordinary differential equations. The details of the derivation are given in the appendix A. 
Explicitly we find 

tr— ^ J = -8^(1 - a)a ( Ha{y)dy . (22) 

The function Ha{x) is the solution of the following set of ordinary non linear differential equations (we 
omit the subscript a and the dependence on x of the variables for notational convenience) 



H' 

K 

c 
b' 



m , 



-m (6i3i2 + uBi) , 
-m (cBi2 + UB2) , 

— 2 m /32 u csc(a;) sin(a;/2) — c (1 — a) csc(a;) (1 + cos(a;)) , 
— 2 m/3i u csc(a;) sin(a;/2) ~ ba csc{x) (1 + cos(a;)) , 

771 1 
-— sec(a;/2) (6/32 + c/3i) + - u tan {x/2) , 



(23) 

(24) 
(25) 
(26) 
(27) 

(28) 



where Bi, B2, B12, (32 are functions of x given in terms of H, Xi, X2, c, b, and u by the following set 
of algebraic equations 



cos(a:/2) 
87ra(l — a) 

sin(a;/2) 
87ra(l — a) 









sin(2;/2) H - m {jSi X2 + 132 Xi) + 2m cos(a;/2) u B12 
- cos(a;/2) H ~ m tan(a;/2) {P1X2 + /32^i 



(29) 



TO 8m{x/2){bB2 + cBi) , (30) 
- f3iB2) + (1 - 2a) cos(a;/2)Bi2 , (31) 



-TOsin(a;/2)(cXi -6X2) + mtan(a;/2)(/32Bi 
-4a(a - 1) - to2(4 - 8/3i/32 +bc + 3u^) 

-4cos(a;) (a(a - 1) + m^(w^ + 1)) + cos(2a:)(&c- u^) , 
(2a- l)?iCOs(a;/2) + mtan(x/2)(/3ic- 6/32) . 



The boundary conditions at x 
H{tt 

X2iTT 

u{tt 

b{TT 

c(7r 



TT are 



= 0, 



r(-a) (cosh (^) Im [V> + g + f )] - f sinh (^)) 
22"/i (cos (2a7r) + cosh(7r^)) r(l - 

Xi{it) I a^(l-a) , 



|r(i 



2 1 



0, 

2i-2'^a(l 



^)|r(i 



2 I 



TOr2(l + a) 

fe(7r) I a^(l-a) , 



(32) 
(33) 

(34) 
(35) 

(36) 
(37) 

(38) 
(39) 



where = v4m^— T and if) is the digamma function. The meaning of the extra variables -Bi, B2, B12, 
Xi, X2, u, b, c, (3x and (32 is the same as in [Tl] and is given in appendix A. The trace in ([^ is regularized 
such that it vanishes when x — tt, where there is no vertex point and no logarithmic term is present in the 
entropies. 

In [20j the partition function for a Dirac fermion on the Poincare disk with boundary conditions 
analogous to psp imposed on a geodesic segment has been written in terms of a solution of the Painleve 
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Figure 1: The solid curves are the functions S2{x) (top) and s^ix) (bottom). The points are obtained by 
numerical simulations in a lattice, and correspond to, from top to bottom, s = si, S2 and S3 evaluated for 
the angles x = 7r/4, tt/2 and 3/4 tt. 

VI differential equation. Though we were not able to find an explicit relation, it is likely that our results 
for H{x) could have also an expression in terms of solutions of these type of equations. 

2.2 The coefRcient of the logarithmic term 

Gathering aU the resuhs together, using eqs. jS]), dH) with A = M, HH), and we arrive at 

the following result for a real scalar 

, , '^8k(n-k) , ,2 r , / 

fr[ " (n - 1) Ji/2 Jx 

where we have made explicit the m dependence of H . This equation, together with (|23ll38p gives our final 
expression for the coefficient of the logarithmic term in Sn for a real scalar field (half the complex scalar 
result), which is in position to be evaluated by solving the ordinary differential equations numerically. 

The functions Sn{x) satisfy Sn{x) = Sni^ir — x) , which is a consequence of the symmetry in the entropies 
S„{V) = Sn{—V) due to the purity of the vacuum state. In the figure (1) we have plotted S2{x) and 53(2:) 
for X € [0, tt]. The values of si — s, S2 and S3 for x = 7r/4, tt/2 and 3/47r obtained by lattice simulations 
are also plotted. They show a perfect accord (less than one percent error) with the analytical results. 
These particular values of the angle are the ones for which the coefficient can be calculated in absolute 
terms with very small error on a square lattice of limited size (in the present case it was 100 x 200 points). 
The numerical methods consist of evaluating the entropy for a massless real scalar (see Appendix B and 
[l4l [21] ) for a given shape (square, triangle, etc.) and different overall size A, and then fitting the result 
as = Co + Ci A + C_i A^^ + C_2 — s„ log(A). It is also possible to evaluate very accurately s„ for 
specific combinations of angles using rectangular triangles. In this way we have computed in the lattice 
Sniu) + Sra(7r/2 — y) with y — arctan(p/(j), where p and q are small integers. We also obtain in this case a 
perfect accord with the analytical results. We have also checked that s„ does not depend on the orientation 
of the polygon with respect to the lattice simmetry axes. 

3 The small angle limit and the two dimensional entropic c- 
functions 

From (|40p and the differential equations for H, it is possible to derive expansions for s„(a;) and s{x) for 
small X or t: — x. In particular, in the small x limit the Green function on the cut sphere is related to the 
corresponding one for the flat space problem with boundary conditions imposed on an interval. Taking into 
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Figure 2: The dashed sets A and B arc anti-starshaped, in the sense that any ray from the origin has no 
intersection with them or the intersection is a half hne which does not contain the origin. A contraction 
of the space transforms them into A' and B' . In the hmit of an infinite contraction they are mapped to 
two angular sectors with common vertex. 



account the results of [14], which deal with this later problem, we can show directly from the differential 
equations that for small x 

six) ^ ^-^^^^ , (41) 
nx 

where cs{t) is the one dimensional entropic c- function for a free scalar field [51[T31[I1]. The c- function for 
a theory in 1 + 1 dimensions is defined as 

c(r) = , (42) 

where S{r) is the entanglement entropy corresponding to an interval of length r in 1 + 1 dimensions. For 
free fields c(r) = c{t) where t — mr, and m is the field mass. 

We can obtain the formula (HT]) (and the analog one which is valid for free fermions) with a less technical 
derivation. This also sheds light on the origin and the necessity of the logarithmic term. In a previous 
paper we have shown that the entropy corresponding to a spatial rectangle in 2 + 1 dimensions with a short 
side L and long side i?, R/L ^ 1, contains the universal term (included in the finite term 5*0 in eq.lU])) 

m 

5univ ^ -kR/L , (43) 

where fc is a dimensionless function of L and the renormalized parameters of the theory. For a free massless 
theory it is 

1 

k^- dt{nsCs{t)+nFCF{t)), (44) 
Jo 

where cg and cp are the one dimensional entropic c-functions for a real scalar and a majorana fermion, 
and ns and np are the multiplicity of the boson and fermion degree of freedom. This gives k ~ 0.039 for a 
real scalar and k ~ 0.072 for a 2 + 1 dimensional Dirac fermion. Now, we may approximately decompose a 
plane angular sector with small angle a; as a union of thin and long rectangles with bigger size as we move 
further from the angle vertex. Thus, there is a term in the entropy for the angular sector which is due to 
the sum of the terms (j43|) for the rectangles (the term (|43|) is extensive in the direction in which lie the 
different rectangles). In the limit of an infinite partition this gives the desired result 



4 Final remarks 

We have shown that there is a logarithmic divergent term with universal coefficient for the entanglement 
entropy corresponding to spatial polygonal sets in 2+1 dimensions. We have found the analytic expression 
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for the coefficient in the alpha entropies for integer a — n and a free scalar field, and the small angle limit 
for the corresponding coefficient on the geometric entropy for general free fields. 

Interestingly, since this paper first appeared as a preprint, a vertex induced logarithmically divergent 
term in three dimensions has also been found in the conjectured geometrical form of the entanglement 
entropy for some conformal field theories arising in the AdS/CFT duality context [H1H2]- The overall form 
of the logarithmic coefficient as a function of the angle is similar to our exact results, and the authors have 
also been able to show the relation between k and s (eqs. (|43p and (US])) for their entropy functions, which 
presumably correspond to some interacting theory. 

Technically, the calculation in this paper amounts to the one of the conformal anomaly for a free 
scalar on a three-dimensional manifold with a conical singularity located on the boundary of a plane 
angular sector. This is a non-trivial result, since there is no known general method applicable to this 
case involving conical singularity surfaces which are themselves non-smooth. We have shown that this 
particular problem can be mapped to the one of the calculation of the trace of the Green function for a 
massive scalar field on a two dimensional sphere where boundary conditions are imposed on a segment of a 
great circle. We have found this trace analytically by a method introduced in [HI [19] which allows one to 
exploit the rotational symmetries of the sphere even in the presence of the symmetry breaking boundary 
conditions. 

Our results for the Green function on the cut sphere may find different uses beyond the present one. 
One possible application is the problem of scattering of waves by a plane angular sector in three dimensions 
[16j . Remarkably, it also gives the exact entropy functions for a spatial segment in 1 + 1 dimensional de 
Sitter space. This space is equivalent to the surface 

r^-t^= i?2 (46) 

in 2 4- 1 Minkowski space. The euclidean de Sitter space then corresponds to a sphere with radius R. The 
formulas corresponding to p2| and ([T3l) in this case take us again to the Green fimction on the cut sphere. 
The result can then be readily written off, 

^de Sitter (^/^) ^ _ g 8nkin- k) r ^ ^ ^^^g ^ (47) 

^ n^{n-l) J^MR)^ Jr/R " 

where r is the physical size of the geodesic interval, and M is the field mass. The additive constant is 
logarithmically divergent with the cutoff, in order to keep Sn positive for small values of r. 

The formulas ([4T|) . ([43]) and ([45]) show a remarkable relation between universal terms in the entropies for 
different dimensions. The function c(r) = rdS(r)/dr, with S(r) the one dimensional entropy corresponding 
to an interval of length r, plays the role of the Zamolodchikov's c-function in the entanglement entropy 
c-theorem [5] . This theorem states that there is a universal dimensionless quantity in two dimensions which 
is decreasing under scaling and has a well defined value at the fixed points (proportional to the Virasoro 
central charge) [23] . The c-function introduced by Zamolodchikov is constructed from a correlator of stress 
tensor traces, and we have shown that it can be taken to be c(r) as well. There has been a substantial 
effort to extend the c-theorem to higher dimensions, but a definitive result in this direction is still missing 
[24]. At fixed points the function c(r) is a constant given by the coefficient of the term proportional to 
log(e) in the two dimensional entropy. Then, one would be tempted to speculate that a running function 
which takes the value s{x) at fixed points could be a good candidate for a c-function in three dimensions. 

On the other hand, the logarithmic term in the entropy induced by the vertices is, remarkably, the only 
obstacle in the following simple argument attempting to prove the c-theorem in any dimensions. Consider 
the mutual information 

IiA,B)^S{A) + S{B)^S{AuB) (48) 

between two non intersecting sets A and B. This is dimensionless and universal (all boundary terms get 
subtracted), and it is also positive and increasing with the size of each set separately [5]. Thus, if we take A 
and B as anti-starshaped (see figure (2)) I{A, B) is decreasing under dilatations. However, it fails to have 
a well-defined value at the ultraviolet fixed point, since in that limit A and B go to angular sectors with 
a common vertex, and I{A, B) has a logarithmic divergence due to the mismatch of the vertex induced 
terms in (|48]) . 
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5 Appendix A: Green function on the cut sphere 

In this appendix we derive the set of non linear differential equations (|23ll38p which give the trace of the 
Green function on a sphere with a cut of angle x (see figure 3), where the boundary conditions (jlSp are 
imposed. We will follow here the same steps as in |14j where we have solved the analogous problem on 
the plane with boundary conditions imposed on a finite interval. We use a generalization of a method 
originally developed in [19]. 

The Green function G{z,z') = (—An + m?)~\,, where z and z' are points on the sphere, is uniquely 
defined by the following three requirements: 

a.- It satisfies the Helmholtz equation on the sphere 

(-Ajj + m2) G(z, z') - 6{z - z') , (49) 

where, in polar coordinates, 

19 d 1 92 

sm B oO oO sm 6 oip-^ 



b.- The boundary condition is (it also holds for the Green function derivatives) 

lim Giiir + e, ip),z') = e'^-" lim G{{n - e, ip),z') for ip e ^i, t/Ja] ■ (51) 

c— ►0+ e— >0+ 

Here the points with polar coordinates (7r/2,(p), tp G [ipi,ip2], form the cut location. We need to specify 
the cut through its two cndpoints with angle ipi and if 2- At the end of the calculation we can set as in 
Section 2 ipi = —(p2 = x/2. From now on we also choose a G [0, 1]. 

c- G(z,z') is bounded everywhere (including the cut) except at z = z' . 

Several properties of the Green function can be derived from its definition. First it is Hermitian 

G{z,z') = G{z',z)* , (52) 

and the two reflection symmetries of the problem give 

G{z,z') = G{Tz,Tz')*, (53) 
G{z,z') ^ G{Rz,Rz'), (54) 

where R and T are the reflexion operations given by 

Ri9,cp) = i9,ip,+ip2-v), (55) 
Ti9,ip) = {7r-e,ip). (56) 

In what follows we study the structure of singularities of some functions related to G{z, z') and use 
repeatedly the fact that a nonsingular solution of the homogeneous Helmholtz equation must be zero due 
to the uniqueness theorem. 

According to the boundary conditions, near the end points of the cut the Green function must have 
branch cut singularities. The requirement that the function must remain bounded on the cut and the 
equation ([49]) imply that the leading terms of G{z, z') for (6*, 93) near (7r/2, (pi) (and fixed z') have to be of 
the form 

G(z, z') ^ \{f - ^1) - i{B - V2)]"5i(z') + \W - ^1) + *(e - ^l2)\^-''S2{z') , (57) 

for some functions S\i^z) and S2{z). We have written explicitly only the terms with powers of \(ip — Lp\) — 
i{9 — 7r/2)] with exponent smaller than one. These are the ones of interest in the sequel. Note that the 
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Figure 3: The plane angular sector of angle x = ipi — ip2 and the sphere with a cut of angle x 



contributions at this order must be analytic or anti-analytic in z in order to cancel the Laplacian term in 

dm). 

The most singular contributions to d^^G{z, z') for z (7r/2, ipi) follow from the derivative of ([57)1 

5^,G(z, z') ^ -a[((^ - ^i) - i{e - ^/2)f-^Si{z') + (a - - ^i) - i{e - V2)]~'^52(z') . (58) 

The function d^-^G{z, z') satisfies the homogeneous Helmholtz equation and the boundary conditions, and 
it is not singular at z — > z'. It has only one singularity located at ipi whose expression is given by (|58p . 
Thus, an adequate linear combination of this function for different values of z' must be a nonsingular 
solution of the Helmholtz equation, and therefore identically zero. Following the same argument as in the 
flat case rRT, this leads to the fundamental relation 



d^,Giz,z') = Siz^ASiz'), 



where we are using vectorial notation 



S{z) 



Si{z) 
S2{z) 



(59) 
(60) 



and A is an hermitian matrix. 

The function S{z) satisfies the homogeneous Helmholtz equation with the same boundary conditions 
as G{z' , z) with z' fixed, but it is unbounded around z = {tt/2, tfi). In fact, it follows from ([58]) and (|59p 
that it has singular terms proportional to [{(p—ipi) — i(0 — 7r/2)]^° and [{ip — ipi) + i{9 — tt /2)]°-''^ . However, 
these disappear if we derive with respect to ip2- On the other hand S{z) behaves around (ixjl^ip^) as an 
ordinary wave, that is, it vanishes proportionally to \(lp — ip^) — i(G — 7r/2)]'' and \(lp — ip^) + «(^^ ~ '''/2)]"'^~"- 
Then d^^S{z) has singular terms around (j^jl^Lp-i) which are proportional to \(lp — (^2) — — 7''/2)]''~^ 
and \{ip — tpi) + i[d — 7r/2)]~°. This structure of singularities and the relation ((54| gives place to 



d^^S{z) = ^S{Rz), 
d^,S{Rz) = --iS{z), 

where the matrix 7 is a function of ipi — p>2- 

Using ([55)1 . and to calculate dip^d^^G{z, z') — d^^d^^G{z, z') = we get 

7^ A'^A^^ , 



(61) 
(62) 



(63) 



and that the matrix A must be a constant, ^A/^^p^ — dA/d^^ — 0. Thus, it can be evaluated by the 
knowledge of the solutions for a; = vr or using the flat space limit a; ^ [H] . It is given by 



A — — 47r(l — a)a ui , 



where cti is the Pauli matrix. 

The equation ([SS)) leads to the behavior 



S{z) 



1 

47r 



(64) 



(65) 
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for z in the vicinity of (7r/2, Lp\). Using we get a similar expression for S(Rz) 

^ ' 47ra(l-a) V V2)] / ^ ' 

5.1 Equations for S[^z) from rotation symmetries 

In order to use eq. ([55)) to compute the trace of the Green function we need more information on S{z). With 
this aim we exploit the symmetries of the problem as follows. We consider two rotations around the axes 
lying on the equator which cut the sphere at the points (Q ~ — ^£i±^) and {9 — t: /2,ip — + 

The associated differential operators are given respectively by 

Li = — cot 9 sm{(p — (p'^)dtp + cos{(p — (p'^ jdO , (67) 
L2 = — cot 9 cos{(p — (p^)dip — sin{(p — (p^)d9 . (68) 

We use the notation (p^ — '^^^'^^ . From ([57)1 and ([S5)l we check that they satisfy the algebra of angular 
momentum operators 

= d^, (69) 
[Li,d^] = -L2, (70) 
[L2,d^] = Li, (71) 
Ll + Ll + dl = An, (72) 



and the following relations 



dLi 

dipi 

dL2 1 



2^- 

They commute with the Laplacian operator and thus the functions {Li + L'^ G{z, z') with i = 1, 2 are 
non singular at z = z' and satisfy the homogeneous Helmholtz equation and the boundary conditions. 
However, due to the LiG{z, z') terms, they are singular at (7r/2, tpi) and (7r/2, 992) as a function of z. Near 
(7r/2, ipi) we have from ([57 )) . ([65)) . ([67)) and ([68| 



LiG(z, z') ^ cos^~ [ia[{(p ~ (p^) - i{0 - Si{z') - i{l - a)[{^ - (^1) + i{0 - tt/2)]-''S2(z')] 

- -icos{ip-)S{z')Aa3S{Tz) , (75) 
L2G(^,^') ism{^~)S{z')Acr3S{Tz). (76) 

This gives from the uniqueness of the solution of the homogeneous Helmholtz equation for non-singular 
functions 

LiG{z,z') + L[G(z,z') = -icos{ip-)[S\z)a3AS{z') + SHRz)a3AS{Rz')] , (77) 
L2G{z,z') + L'^G(z,z') = -ism{(p-)[S^z)a3AS{z')^ S''{Rz)cr3AS(Rz')] . (78) 

From ([65)) we see that the most divergent terms of LiS^{z)A — i cos{(p^)dipS^ {z)(73A exactly cancel. 
But this combination should also have a contribution to order [{ip — (^2) — — ''■/2)]"^^ and [((/s — (^1) + 
i{9 — 7r/2)]~° for z in the vicinity of (7r/2, Lp\). Thus we have 

LiS\z)A - ico^{ip-)d^S\z)a3A = ~S\z)A^,, (79) 

where ^ is a matrix which depends ow x ~ Lp\ — (^2- A similar argument hold for L2S'^{z)A + i sin(<y9^) 
d^pS'^ {z)g3A giving place to 

L2S\z)A^is\x^{ip-)dipS\z)a3A^S\z)Aix, (80) 
with ^1 depending on a; = (/Si — (^2- 
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Taking the derivative with respect to ipi of ([77)1 and ([75)) . and extracting the coefScients of the divergent 
terms, we get 

LiS'(z) = CS{z)-ism{ip-)a3Siz)-icoa^-[{j,a3}S{Rz)-a3d^S{z)] , (81) 
L2S{z) = -^iS{z)-icos{ip-)a3S{z)+ism{^-)[[a3,j]S{Rz) + a3d^S{z)] . (82) 

We can also write the reflected equation for ((8T|l and ((82|) 

LiS{Rz) = CS'(i?^) - isin(p-cr3S'(i?z) - icos<y9~ [{7, erg} S'(z) - cT3,d^S{Rz)] , (83) 
L2S{Rz) = Ci'S'(-R^) + j cos((^")cr3S'(i?z) - i sin((^-) [[0-3, 7] S{Rz) - a^d^Siz)] . (84) 

Subtracting the hermitian conjugate of the equation (|5T|) and ([55]) and using {A, (T3} = 0, {7, 173}^ A + 
^{li ""s} = Oj obtain 

£}A + A£, = -isintp^a^A. (85) 
A similar relation for ^1 follows by subtracting the hermitian conjugate of ([5^ and (|84p . 

+ = icos(p~CT3-4. (86) 



5.2 Parametrization 

The expansion of S{z) can be extended to the following order around the singular point (7r/2,<y9i) intro- 
ducing a real matrix TV 

<^(z\ ^ ^ ( -aiiv-Vi)- ^{e - 7r/2)]- \ / [(^ _ _ ,(0 _ ^/2)]i- \ 



47r V (T3^[(^ - ^1) + - ^1'^)] ) \ [(^ - ^1) + ^{d - TT/2)f 



This general expansion for S inserted in (I79p and (|80p allows us to solve for some components of N and 
obtain a relation between ^ and ^1 

^1 = tan((^~)^ + igsec((^~) , (88) 

with 

■-'J. (89) 

It also gives additional relations which, together with the algebraic equations ([63]), (|85p . give the general 
form of 7 and ^ in terms of the following parametrization 

m/u fo\ _ / -i(a- l)sin((^-) imfii \ 

^ 2\c u ) ' ^ -iml32 -msin((^-) ^ ' ^ 

where u, 6, c, /3i, and (32 are real functions of m and 

Taking derivatives respect to 9 and if of (|8ip and ([5^ and combining them to reconstruct the Helmholtz 
equation for S{z), and using the expansions of S{z) and S{Rz), we get from the series around the singular 
points the differential equations ((26)) . ([27| . ((28| for 6, c and u, and the algebraic equations ([32)1 and (|33p 
which determine /?i and /32- 

5.3 Integrated quantities 

The equation ([59]) gives for the trace of the Green function 



^ J dfl G{z, z) = J dfl S-U S = ~8na {1 - a) H , 



where 

H{x) = / dnS*iiz)S2{z), (91) 

and dfl = d9dipsm{9). 
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To find this quantity, we use the information obtained in 5.1 and 5.2. We first define the following 
auxiliary integrals 

(92) 
(93) 
(94) 
(95) 
(96) 

These are convergent. They are also real since the relation ([54)) implies that S*{z) ~ Si{Tz). 

The equations for these quantities are obtained basically combining conveniently the components of 
([5T|) and ([5^ multiplied by the components of S{z) and S{Rz), and integrating on the sphere. In this way 
the combination {Li Si{z)) 5*1(2:) + >S'i(2;) {Li 6*2(2))* integrated on the sphere gives place to the equation 
Similarly, from the combination (L2'5*i(z)) 5|(z) + Si{z) {L2S2{z))* we get equation ([5D|) and from 
{L2S2{Rz)*) Si{z) + S2{RzY (L25i(z))* we get equation 

Differential equations for the integrated variables follow by taking the derivative of ipij) , ([55)1 and (|96p 
with respect to ip2- These correspond to (|23p . ([M)) and ([^5]) of Section 2.1. 

In the flat space analog of the present calculation it is possible to give a closed algebraic expression of 
H in terms of u [14j . This simplification seems to be lost in the sphere. 





/ <70 '?*('7'i Si 


iJ2\X) — 


/ WO Q*( r>\Q I 


Bl2{x) = 


J dns;{z)Si{Rz), 


Xx{x) = 


J dnsi{z)Siiz), 


X2{X) = 


J dns;{z)S2{z). 



5.4 Boundary conditions at x = ir 



The homogeneous Helmholtz equation on the cut sphere for a function f{z) can be solved exactly by 
separation of variables when x = it. The solution has the general form 



(97) 



where we have put now the end points of the cut at the poles 6 — and 9 = n oi the polar coordinates. 
Here /„ is given by 

/„(0) = c^prix^^ iw) + C2Q||:,^,^) H , (98) 

with Q and P the standard Legendre functions, w = cos 9 and /i = y/Am'^ — 1. Evaluating the limits w ^ 1 
and w —1 and comparing with the expansion (j65p for S'(z) we get for x = n 

2-a/2 ^ 



Si 
S2 



47ra 2 2 ^ cos(a7r — csc(a7r)r(a) 

2(a-l)/2 I 

47r(l - a) 2^T^ cos(a7r + csc(a7r)r(l - a) 



-cot(a^)P|(_^^^^)(^)-Q|(_^^^^)(u;) (99) 



-cot((l-a)7r)P|-_\^^^^(u;) 



Jip{a-1) 



(100) 



Using these expressions, the eqs. (IFTt and ([62]), and the expansion of S{Rz) we find u(7r), 6(7r) and c(7r) 
given in dSZl), dSll) and dSH). The results dSl, (ES]) and dS!]) for iJ(7r), Xi(7r) and X2{it) follow directly 
from ([9T|). (l95|) and ((96|) . and the explicit form of 5*i and 5*2, where we have used some formulae for the 
integrals of a product of two Legendre functions given in [25] . 



6 Appendix B: numerical entropy in a two dimensional lattice 

We use the method presented in [21j to give an expression for pv in terms of correlators for free bosonic 
discrete systems. Take a free Hamiltonian for bosonic degrees of freedom with the form 
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where (pp and T:p obey the canonical commutation relations [(j)p^-Kpi] — i6pp, M is a Hermitian positive 
definite matrix, and the sums are over the lattice sites r. The vacuum (ground state) correlators are given 

by 

Xpp, = = i(M-^),>-, , (102) 

PPP, = {TipTip,) ='^{A&)pp, . (103) 

Let XYp, and Ppp, be the correlator matrices restricted to the region V, that is r, (E V . The entropies 
can be evaluated as 

S ^ ((^'^ + 1/2) log(^^fe + 1/2) - {uk - 1/2) log(i.fe - 1/2)) , (104) 
fc 

Sn = V log ((i^fc + 1/2)" - - 1/2)") , (105) 

n — 1 ^ — ' 
fc" 

where Vk are the eigenvalues of VX^.P^. Thus, to compute the entropy numerically we need to diagonalize 
the matrix X^.P^. 

We took the lattice Hamiltonian for a real massless scalar in three dimensions as 

1 °° 

= 2 X! ('^"™ ^4>n+l,7yi - (t'n,mf + ('/'n,m+l " </>ri,m)^) ■ (106) 
n,m— — oo 

We have set the lattice spacing to one. The correlators (|102p and (jlOSp are 



^ r r cos(a:(^~l))cos(j/(J^-l)) 
87r2 ^ ^ v/2(l - cos(x)) + 2(1 - cos(y)) ' 



i'oo'/'y) = I / dxdy ^— ^ =, (107) 



{^^oo^^^J) ^ / / da^dy cos(a;(i - 1)) cos(y(j - 1)) ^2(1 - cos(a;)) + 2(1 - cos(2;)) . (108) 



TT /•TT 



— TT — TT 



It is relevant to the accuracy of the entropy calculation to evaluate the correlators with enough precision. 
This can be very time-consuming. We have found it is much faster to evaluate one of the two integrals in 
the correlators analytically (in terms of polynomials times elliptic functions) using a program for analytic 
mathematical manipulations, and then doing the last integral numerically. 
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